# power for t corr

r=0.3; n=8;

d=sqrt(r^2*(n)/(1-r^2))

tc = tinv(0.95,n-2)

power = 1-nctcdf(tc,n-2,d)